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Magic wavelengths for laser trapping of boson isotopes of alkaline-earth Sr, Ca and Mg atoms are investigated 
while considering terahertz clock transitions between the 3 Pq, 3 Pi, 3 Pi metastable triplet states. Our calculation 
shows that magic wavelengths of trapping laser do exist. This result is important because those metastable states 
have already been used to realize accurate clocks in the terahertz frequency domain. Detailed discussions for 
magic wavelength for terahertz clock transitions are given in this paper. 

PACS numbers: 0.6.20-f; 0.6.30-Ft; 32.80.Qk; 32.30.Bv. 



I. INTRODUCTION 
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Frequency standards have achieved an unprecedented suc- 
cess in experimental demonstrating accuracies of 4 x 10~ 16 
with a cesium microwave fountain clock Ut] and 1.9 x 10 -17 
with an ion optical clock J2, Ht] . For optical frequency stan- 
dards based on neutral atoms, in order to effectively in- 
crease the interrogation time, Katori proposed to utilize op- 
tical lattice trap formed with a magic wavelength trapping 
laser 0, @], This clever technique greatly enhanced estab- 
lished high-accuracy optical frequency standard with neutral 
Sr atom to an accuracy of 1 x 10" 16 
clock schemes based on Ca U, Yb \\ 
with magic wavelength lasers have been proposed. 

Optical trap with a far off -resonant laser is a very useful tool 
for the confinement of cold atoms. Nevertheless, for the preci- 
sion of clock transitions in frequency standards, light shift due 
to trapping laser has to be avoided. Thus the wavelength of 
the trapping laser should be tuned to a region where the light 
shift for the clock transition is eliminated, that means the light 
shifts of the two clock transition states cancel each other. The 
wavelength A is called magic wavelength BHH]. Recently, ce- 
sium primary frequency standard with atoms trapped in an op- 
tical lattice with a magic wavelength was suggested lfl2l[l3ll . 
and possible magic wavelengths for clock transitions in alu- 
minium and gallium atoms were also calculated Hill . 

In contrast to the above mentioned magic wavelengths for 
optical clock transitions and microwave clock transitions, here 
we investigate magic wavelengths for terahertz clock tran- 
sitions. Absolute frequency standards in the terahertz do- 
main with fine structure transition lines of the Mg and Ca 
metastable triplet states were first proposed in 1972 by Stru- 
mia fl5ll . After more than twenty years of continuing im- 
provement, a frequency standard based on the 3 J°i - 3 Pq Mg 
transition and thermal atoms in a beam has reached an un- 
certainty of 1 x 10~ 12 flril [nil . However, these potential 
terahertz transitions for high-resolution clock references have 
never been experimentally investigated with laser cooled or 
laser trapped atoms. 

In this paper, we present our most recent calculation 
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FIG. 1 : (Color online) Simplified diagram of the lowest energy levels 
for alkaline-earth atoms and some possible laser couplings. 



of trapping laser magic wavelengths for Sr, Ca and Mg 
atoms, considering different possible clock transitions be- 
tween metastable triplet states 3 P. Accurate terahertz clocks 
could then be built based on such atoms which are cooled and 
trapped in an optical lattice. 



II. THEORETICAL DESCRIPTION 

For alkaline-earth atoms, two valence electrons result in 
two series of atomic energy levels as the electron spins can 
be parallel (triplet states) or anti parallel (singlet state). The 
energy diagram can be simplified as shown in Fig. 1. The 
ground state is So, and the lowest excited states nsnp are 
l P\ and 3 Pj which can be divided into three fine structure 
sublevels 3 P2 3 P\ 3 Po- For the 3 Pj state, transitions to higher 
states can be divided into three groups : 3 Pj - 3 S \ 3 Pj - 3 Pj 
and 3 Pj - 3 Dj. 

By second-order perturbation theory, the energy shift 
Uj(a>, p, ntt) of atomic state |/) with energy and Zeeman 
sublevel m,, which is induced by a trapping laser field with 
frequency v-to/ln, polarization p, and irradiance intensity / 
can be expressed as Uj(a>, p, mj)=-ai(u>, p, mi)I/2eoc with the 
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induced polarizability or,-. 

The polarizability can be calculated by summing up the 
contributions from all dipole interactions between the fine 
structure state |/) and \k) with the Einstein coefficient A jki 
(spontaneous emission rate for Ek>Ef), Zeeman sublevels 
OTi, rtik and transition frequency vjki=cojki 
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Here e is the electron charge, ha> jki is the energy difference be- 
tween fine structure states \k) and |/), /3 denotes other quantum 
numbers of the state, and (fikJk \\D\\/3iJi) is the dipole reduced 
matrix element. The expression in large parentheses in Eq.(l) 
denotes a 37 symbol which describes the selection rules and 
relative strength of the transition depending on the involved 
angular momenta J, the projection m, and the polarization p. 



If we know aijki and A/& in Eq.(l), we can get the polar- 
izability cc,. However, typically the literature gives the total 
transition rate A T from a given excited state to the fine struc- 
ture manifold states below. So we need establish the relation 



TABLE I: Sr element: Transition Wavenumbers (WN)(cm ') corresponding to (Ojki, Einstein Coefficients for fine structure states Aj(xl0 6 s ') 
and Total A T (xl0 6 s~ [ ) for fine structure states manifold, Correction Factors f . The Wavenumber data originate from (20ll . 
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TABLE II: Ca element: Transition wavenumbers (WN)(cm 4 ) corresponding toa> Jkj , Einstein Coefficients for fine structure states A j(xl0 6 s~ l ) 
and Total A T (xl0 6 s^ 1 ) for fine structure states manifold, Correction Factors f . Besides the updated values listed in Ref. 13311 for A T , the others 
originate from l33l . 
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between Ajki and A j^i- We know Ajm can be expressed as: Here 
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Here fiurki is the energy difference between two fine structure 
manifold states \k) and \i). Using the formula: 



(J3 k J k WDWPiJt) = (-l) Lt+St+Ji+1 S StSi y/(2Jk + 1)(2// + 1) 
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and combining Eq.(2) and (3), we can get 
Ajki = A Tki x ^{(jj ki )R ki 



(5) 
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is the energy dependent correction lfl9h . reflecting the alter- 
ation on the transition rate due to the effects such as the orbit- 
spin interaction and the spin-spin interaction which causes the 
fine structure splitting. And 



Rki = (2Z* + l)(2/j + 1) X 



Jk 1 Ji 
Li Si L k 



(7) 



gives the fraction of the coupling strength between an excited 
state and a lower state | i). Since the total transition rate Ajki 
is usually available in the literature, this geometric ratio tells 
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us how to scale the interaction for a particular fine structure 
state of interest. 

To calculate the wavelength dependent polarizability, we 
combine Eq.(l) with Eq.(5), and use the known transition fre- 
quencies and spontaneous emission rates in the literature. This 
light polarizability is very sensitive to the Einstein coefficient. 
However, theoretical and experimental values of magic wave- 
length for the optical clock transition obtained in the past can 
be used to confirm our calculation. 

In this paper, we use this method to calculate the light shift 
for the terahertz clock transition from 3 Po to 3 Pi,m = lev- 
els, and from 3 P\,m = to 3 P2,m = levels for boson iso- 
topes with the nuclear spin 7 = 0. After calculating magic 
wavelengths for Sr and Ca optical clock transitions and com- 
paring them to the experimental values, we calculate the po- 
larizability of terahertz transition with data collection mainly 
from Ref . J2C " 



HI. CALCULATION OF MAGIC WAVELENGTH 
A. Strontium 

Using the method above, for Strontium, we first calculate 
the magic wavelengths of two optical lattice clock transitions 
with the datas listed in Table I and compare the results with ex- 
perimental values. Then, we calculate the crossing points for 
terahertz clock transitions where the difference of polarizabil- 
ity is zero. Table I shows Transition Wavenumbers, Einstein 
Coefficients and Correction Factors for the 5s 2 l So, 5s5p 3 Pq, 
5s5p 3 Pi and 5s5p 3 P2 states for Sr element. For Einstein 
Coefficient Atm, first we choose the available updated exper- 
imental values in Re f. 11261 128113 ill , then we use updated the- 
oretical data in Ref. ||27| 1 3 Oil , and for the rest we mainly use 
theoretical values in Ref. y5[] . 

According to our calculation, the crossing point for the 
1 So to 3 Po transition occurs at 813.1 nm, while the crossing 
point for the So to 3 P\(mj = +1) transition with linear po- 
larized light takes place at 915.4 nm. Both of those results 
are in a gree ment with the experimental values of 813.428(1) 
nm h21| - |24|1 and 914(1) nm 112711 . This confirms our calculation 
procedure. 

Fig. 2 and Fig. 3 display the wavelength dependence of the 
atomic polarizability difference Aa for Sr with trapping laser 
wavelength around 400 nm and 1650 nm, respectively. The 
result is scaled by a factor of l/(4nsaal) and the polarizability 
is given in atomic unit. In Fig. 2, for linear polarized light, 
Aa between 3 Pi and 3 Po and Aa between 3 P2 and 3 Pi are 
given in solid and dash dotted lines, respectively. In Fig. 3, Aa 
between 3 Pi and 3 Po for linear polarized light and Aa between 
3 P2 and 3 Pi for circular polarized light are presented. The 
cross markers are the crossing points where Aa is zero. From 
Fig. 2 and 3, we can know that the magic wavelength for 3 Pq 
to 3 Pi,m = with linear polarized light are 381.2 nm, 413.6 
nm, 419.3 nm, 1714 nm and 3336 nm, while for 3 Pi,m = 
to 3 P 2 ,m = are 384.5 nm,441.9 nm and 511.0 nm. On the 
other hand, for m — and circular polarization of light, the 
magic wavelengths for 3 Pq to 3 P\,m = transition are 511.8 
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FIG. 2: (Color online) The wavelength dependence of the difference 
of atomic polarizability for Sr element around 400 nm. 
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FIG. 3: (Color online) The wavelength dependence of the difference 
of atomic polarizability for Sr element around 1650 nm. 



nm and 662.8 nm, while for 3 Pi,m — to 3 P2,m = 0, the 
magic wavelengths are 717.7 nm and 1591 nm. 

B. Calcium 

We calculate the polarizabilities using the data in Table II 
with the same method. Table II shows Transition Wavenum- 
bers, Einstein Coefficients and Correction Factors for the 
4s 2 l So, 4s4p 3 Po, 4s4p 3 P { and 4s4p 3 P 2 states for Ca. For 
Einstein Coefficient, we use the updated theoretical values ac- 
cording to Ref. 13311 . and others are from the data listed in 
Ref. 13211 . In order to check the accuracy of our calculation 
and the data used, we get the magic wavelength 799.2 nm for 
the ^o, m = to 3 P\,m = optical transition with circularly 
polarized trapping light, which agrees well with the experi- 
mental value 800.8( 22 ) nm (H. 

The wavelength dependence of the atomic polarizability 
difference Aa around 350 nm and 1350 nm are shown with 
atomic unit in Fig. 4 and 5, respectively. The crossing points 
where the Aa is zero are marked by cross. The magic wave- 
lengths for linear polarization occur at 1361 nm and 2066 nm 
for clock transition 3 Pq - 3 Pi , while for the transition between 
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TABLE III: Mg element: Transition wavenumbers (WN)(cm"') corresponding to toju, Einstein Coefficients for fine structure states 
Ay(xl0 6 s _1 ) and Total Ar(xl0 6 s _1 ) for fine structure states manifold, Correction Factors £ for Mg element. The Wavenumber and A T data 
originate from l32ll . 
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26045.849 


40.157 


0.9977 


161 


3s3rf 3 D 3 
3s4d 3 D l 




0.000 
28.106 






0.000 
21.040 




26045.867 
32281.157 


160.630 
1.397 


0.9977 
0.9981 


161 
50.4 


32341.930 


1.0038 


32321.871 


1.0019 


3sAcPD 2 




0.000 




32321.830 


37.872 


1.0019 


32281.116 


12.576 


0.9981 


50.4 


3s4cPd 3 
3s5d 3 D, 




0.000 
13.101 






0.000 
9.808 




32281.078 
35057.093 


50.304 
0.652 


0.9981 
0.9983 


50.4 
23.5 


35117.866 


1.0035 


35097.807 


1.0017 


3s5d 3 D 2 




0.000 




35097.784 


17.655 


1.0017 


35057.070 


5.865 


0.9983 


23.5 


3s5d 3 D } 
3s6d 3 Di 




0.000 
6.967 






0.000 
5.217 




35057.040 
36531.700 


23.460 
0.347 


0.9983 
0.9983 


23.5 
12.5 


36592.473 


1.0033 


36572.414 


1.0017 


3s6d i D 1 




0.000 




36572.381 


9.391 


1.0017 


36531.660 


3.120 


0.9983 


12.5 


3s6d } D 3 




0.000 






0.000 




36531.657 


12.479 


0.9983 


12.5 











T 



level P\,m — and P2,m — at 312.2 nm, 316.2 nm, 
325.4 nm, 344.0 nm and 393.4 nm. 

The laser polarization have no effect on the polarizability 
for the ground state (7 = 0) because the ac Stark shift is iden- 
tical with any polarizations. It is also true for 3 Pq state. How- 
ever, the influence of circular polarized laser light is worth 
study for other states. For m = 0, we can obtain the magic 
wavelengths for the 3 Pq to 3 Pi clock transition at 301.0 nm 
and 310.0 nm, while for 3 P\ to 3 P2 one finds 1318 nm and 
2254 nm. 



C. Magnesium 

With the completion of the NIST database, the atomic po- 
larizability of the Mg triplet states in the presence of linear 
and circular polarized light can also be calculated. Table III 
presents Transition Wavenumbers, Einstein Coefficients and 
Correction Factors for the 3s 2 1 Sq, 3s3p 3 Pq, 3s3p 3 P\ and 
3s3p 3 P2 states for Mg element. Using the data presented in 
Table III, the magic wavelengths of 3 Pq to 3 P\,m = transi- 
tion with linear polarization are 335.6 nm and 399.5 nm. The 
magic wavelengths are 308.6 nm, 336.5 nm, 406.1 nm for the 
transition between 3 P\,m = and 3 P2, m = 0. 

For circular polarization of light, the magic wavelengths for 
the transition 3 P to 3 P u m = are 307.7 nm, 336.4 nm, 407.8 
nm. However, we can not find any magic wavelength for cir- 





3 p.- 3 p.. a 



300 310 320 330 340 350 360 370 380 390 400 

wavelength X (nm) 

FIG. 4: (Color online) Wavelength dependence of the difference be- 
tween excited state and ground state atomic polarizability around 350 
nm for Ca element. 



cular laser between level 3 P\,m = and 3 Po, m = 0. 

For Mg atoms, several optical transitions between the en- 
ergy levels of terahertz clock transition states and other levels 
exist, such as 456.5 nm ( 3s 2 l S - 3s3p 3 P ), 383.6 nm ( 
3s3p 3 P - 3s3d 3 D ), 309.6 nm ( 3s3p 3 P - 3s4d 3 D ), 333.2 
nm ( 3s3p 3 P - 3s5s 3 S ) and 517.4 nm ( 3s3p 3 P - 3s4s 3 S 
). Hence, not all the magic wavelengths are good enough for 
clock transition, because the slope of the light shift difference 
with the wavelength is too large (shown in the final table IV). 
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FIG. 5: (Color online) Wavelength dependence of the difference be- 
tween excited state and ground state atomic polarizability around 
1350 nm for Ca element. 




300 320 340 360 380 400 420 

wavelength X (nm) 



FIG. 6: (Color online) The wavelength dependence of the difference 
of atomic Polarizability around 400 nm for Mg element. 



To some extent, a possible magic wavelength near 400 nm is 
shown in Fig. 6 with the atomic unit. In Fig. 6, Aa for 3 Pi - 3 Po 
transition and 3 Pi - 3 Pi transition with different polarization 
are given. The cross markers reflect the crossing points where 
the atomic polarizability difference is zero. 



IV. DISCUSSIONS AND CONCLUSIONS 

In summary, we have calculated magic wavelengths for ter- 
ahertz clock transitions for alkaline-earth atoms. The calcula- 
tion results are presented in Table IV along with the slopes 
of the difference of polarizabilities at corresponding magic 
wavelengths. Depending on the calculation and the current 
laser development, we recommend 1714nm and 1591nm for 
Sr terahertz clock, 1361nm and 1318nm for Ca terahertz 
clock, 399. 5nm and 407. 8nm for Mg terahertz clock, because 
the difference of polarizabilities have small slopes at these 
magic wavelengths, where we ignore the effect of highly ex- 
cited states and continuum states which can only make little 
contribution to the wavelength dependent polarizabilities at 
terahertz region. 



TABLE IV: Magic wavelengths for terahertz region. LI is the linear 
laser for the 3 Po to 3 Pj clock transition while CI is for the circu- 
lar laser, L2 is the linear laser for the 3 Pi to 3 P 2 clock transition 
while C2 is for the circular laser, k is the slope of the shift difference 
of two level of clock transition levels at the corresponding trapping 
laser wavelength with unit Hz/nm, the sign denotes the direction of 
the change with the shift for the high level minus the lower level. 
The data are given in the reasonable experiment condition with in- 
put power 150 mW focused to a waist of 65 fim as the light intensity 
1.1301 * 10 3 W/cm 2 . 



3p 

0,1,2 


Mg 


Ca 


Sr 


A 
(nm) 


K 

(Hz/nm) 


A 
(nm) 


K 

(Hz/nm) 


A 
(nm) 


K 

(Hz/nm) 


Ll 


335.6 


-1125 


1361 


-3.201 


381.2 


-615.7 


399.5 


-103.5 


2066 


54.94 


413.6 


32.22 










419.3 


-31.4377 










1714 


-5.590 










3336 


9.122 


L2 


308.6 


-26574 


312.2 


8319 


384.5 


1792 


336.5 


1905 


316.2 


5879 


441.9 


5759 


406.1 


220.7 


325.4 


1641 


511 


1697.63 






344.0 


542.9 










393.4 


1787 






CI 


307.7 


-3174 


301.0 


8610 


511.8 


252.8 


336.4 


469.9 


310.0 


-4072 


662.8 


-1068 


407.8 


54.76 










C2 






1318 


-3.365 


717.7 


2692 






2254 


13.39 


1591 


-5.551 



In this paper, we are only focusing on the study of pos- 
sible magic wavelengths of trapping laser for these terahertz 
clock transitions of Sr, Ca and Mg atoms. These terahertz 
clock transitions were first proposed as early as 1972 lfl5ll . 
and recently have been proposed to be applied in active opti- 
cal clock 13411 . These clock transitions of alkaline-earth atoms 
correspond to a 0.6 THz to 11.8 THz frequency region. Af- 
ter the successful developments of microwave fountain fre- 
quency standards, optical clocks with trapped ions and optical 
lattice trapped neutral atoms, it is interesting to study clock 
transitions at terahertz wavelengths. The advantages and dis- 
advantages of terahertz magic atomic clock will be discussed 
elsewhere. The wavelength range studied in this paper (from 
500 /mi to 25 pm) corresponding to THz frequency standards 
fills the gap between microwaves and optical waves. 
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